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摘 要 : 台 特 玛 湖 干 酒 湖 盆 区 风蚀 沙漠 化 快速 发 展 ,发 育 了 以 新 月 形 沙丘 和 风蚀 坑 ( 主 要 为 梭 状 坑 
和 和 碟 状 坑 ) 为 主 的 风沙 地 魏 。 风 蚀 坑 的 发 育 不 仅 受 来 流风 影响 ,还 与 坑内 三 维 流 场 有 关 。 基 于 风 
洞 实验 和 计算 机 流体 力学 (Computational fluid dynamics, CFD ) 数 值 模拟 ,对 风蚀 坑 三 维 流 场 进行 了 


探究 。 结 果 表 明 :(1) 从 入 风 侧 到 出 风 侧 , 模 状 风蚀 坑 和 碟 状 风蚀 坑 底 本 
的 风速 变化 呈 加 速 一 减速 一 加 速 变化 模式 ,上 口气 流 压 


减速 一 加 速 的 变化 模式 ,上 口 所 在 水 平 


的 风速 均 呈 减速 一 加 速 一 


强 也 相应 地 发 生变 化 , 槽 状 风蚀 坑 和 碟 状 风蚀 坑 整 体 呈 低压 一 高 压 一 低压 变化 模式 ,其 中 樟 状 坑 
变化 更 为 明显 。(2) 两 种 形状 风蚀 坑 对 气流 均 具有 一 定 的 吸附 效应 ,使 坑 外 一 定 范围 的 风沙 流 被 吸 
入 坑内 ,聚集 能 量 和 风沙 流 , 加 剧 风 蚀 坑 风蚀 发 育 。(3) 风蚀 坑 加 剧 了 地 表 侵 蚀 , 侵 蚀 程 度 与 风蚀 坑 
尺度 有 关 , 风 蚀 坑 尺度 越 大 ,对 地 表 侵蚀 越 强 。 本 文 直观 展示 了 两 种 典型 风蚀 坑 中 的 三 维 速度 、 压 
力 、 风 沙 流 流向 分 布 ,揭示 了 风蚀 坑 的 风沙 流 吸 附 机 制 ,研究 结果 可 加 深 对 风蚀 坑 形 态 动 力学 的 理 


解 ,也 可 为 内 陆 干 酒 湖 盆 沙漠 化 防治 提供 理论 指导 


关 键 il: 风蚀 坑 ; 风 洞 实验 ; CFD 数值 模拟 ; 风沙 流 


文章 编号 : 1000 - 6060(2022)06 -= 1784 - 11(1784~1794) 


沙漠 化 是 全 球 最 严峻 的 问题 之 一 ,严重 影响 人 
民 的 身体 健康 和 生命 财产 安全 。 风 蚀 坑 是 沙 地 和 
沙漠 化 地 区 常见 的 一 种 风沙 地 貌 类 型 ,是 沙漠 化 的 
一 种 发 生 过 程 和 具体 表现 形式 。 风 蚀 坑 是 荒漠 化 
地 表 风 蚀 较 为 严重 的 表现 形式 ,造成 地 表 景 观 破 
碎 ; 在 其 发 育 过 程 中 不 断 释 放风 成 沙 ,为 其 他 风蚀、 
风 积 地 貌 的 形成 与 发 育 提供 了 条 件 , 并 不 断 推动 地 
表 沙 汉化 进程 中, 表征 了 植被 退化 是 风沙 活动 增强 
的 直接 诱因 ”。 

风蚀 坑内 的 气流 运动 是 风蚀 坑 形 态 塑 造 的 直 
接 动 力 ,其 主要 受 自 身 性 质 以 及 来 流 方 向 的 控制 ”， 
形成 独特 的 流 场 结构 ”。 运 用 风 洞 实验 和 数值 模拟 
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对 风沙 地 貌 形态 动力 学 过 程 探 究 是 一 项 较为 成 熟 
的 技术 。 蔡 东 旭 等 ”通过 风 洞 实验 探究 了 沙丘 与 气 
流 之 间 的 相互 作用 ,发 现 新 月 形 沙 丘 与 风沙 活动 之 
间 存 在 相互 反馈 作用 。Liu 等 “运用 数值 模拟 探究 
了 横向 沙丘 和 金字 塔 沙 丘 的 风沙 流 场 特征 。 
Jungerius 等 "通过 数值 模拟 建立 二 维 风 蚀 坑 模型 ， 
模拟 研究 了 风蚀 坑 的 二 维 流 场 变化 情况 。Smyth 首 
次 应 用 计算 机 流体 力学 (Computational fluid dynam- 
ics,CFD) 软 件 对 风蚀 坑内 外 气流 特征 进行 模拟 ”， 
并 根据 实际 情况 对 实地 风蚀 坑 上 口水 平面 的 气流 
进行 了 数值 模拟 与 实地 实验 相对 比 ” ,实现 了 风蚀 
坑 的 气流 三 维 立体 展示 同时 量化 了 大 风沙 侍 暴 事 
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件 。 同 时 , 现 有 学 者 已 对 海岸 . 沙 质 草 原 风蚀 坑 的 
形态 ”气流 ”人 蚀 积 速率 ”和 演化 过 程 “ ”进行 了 
大 量 研究 。 

受 全 球 气 候 变 暖 影响 ,中 亚 地 区 鞠 漠 化 问题 日 
趋 严峻 ,加 之 人 类 对 水 资源 的 不 合理 开发 利用 ,很 
多 内 陆 湖 泊 逐 渐 萎 缩 退 化 ,其 至 干 润 (如 位 于 新 驻 
塔里木 盆地 东部 的 台 特 玛 湖 )。 台 特 玛 湖上 自 1972 年 
干 润 后 , 干 油 湖 盆 风 蚀 沙 漠 化 快速 发 展 ,当地 生态 
环境 恶化 ,严重 制约 区 域 经 济 社会 可 持续 发 展 。 迄 
今 为 止 ,人 台 特 玛 湖 干 润 湖 盆 已 经 历经 半 个 世纪 的 沙 
漠 化 发 展 , 干 润 湖 倪 上 已 发 育 了 多 种 风沙 地 貌 ,其 
中 风蚀 坑 是 主要 的 风蚀 地 貌 类 型 ,形态 各 异 且 分 布 
广泛 ,与 新 月 形 沙丘 ( 链 )、 灌 从 沙丘 共同 构成 了 风 
蚀 沙 汉化 初期 的 风沙 地 貌 景 观 。 该 区 发 育 的 风蚀 
坑 与 海岸 和 沙 质 草原 的 形态 有 所 不 同 ,普遍 不 发 育 
坑 后 积 沙 体 ,地表 物 质 为 河 湖 相 沉积 , 呈 沙 土 和 黏 
土 互 层 结构 , 沙 层 较为 松散 ,以 极 细 沙 和 粉 沙 为 主 ， 
细 沙 含量 次 之 ” ,夹杂 少量 干枯 的 树 柳 芦苇 等 植 
物 的 根系 。 目 前 ,关于 内 陆 干 油 湖 盆 区 风蚀 坑 的 形 
态 动力 学 过 程 研究 鲜 有 报道 。 

台 特 玛 湖 干 润 湖 倪 区 正在 成 为 库 鲁 克 塔 格 沙 
漠 的 新 扩展 区 , 查 明 该 区 风沙 地 貌 形 成 发 育 规律 是 
科学 治理 的 重要 前 提 。 虽 然 一 些 学 者 运用 CFD 软 
件 已 对 风蚀 坑内 气流 进行 了 模拟 ,但 风蚀 坑 的 三 维 
风沙 流 模 式 认识 仍 不 全 面 "。 本 文通 过 风 洞 实验 并 
辅 以 CFD 数值 模拟 , 探 明 内 陆 干 油 湖 盆地 区 典型 形 
态 风 蚀 坑 的 三 维 流 场 特 征 及 对 地 表 侵 蚀 影 响 ,研究 
SAS HY DAE E A Be PP Hd ,为 沙 
漠 化 土地 治理 提供 科学 基础 。 


1 研究 区 概况 


人 台 特 玛 湖 位 于 库 鲁 克 塔 格 沙漠 与 塔克拉玛干 
沙漠 之 间 ( 图 1) ,是 沙漠 中 的 生命 湖 , 是 若 羌 县 取水 
源 地 ,由 塔里木 河和 车 尔 臣 河 汇流 形成 。20 世 纪 60 
年 代 之 前 ,塔里木 河和 车 尔 臣 河 上 游 水 资源 充足 尚 
未 开发 利用 ,大 量 河 水 汇 人 台 特 玛 湖 ,使 湖面 面积 
一 度 达 到 约 150 km 。 但 随 着 上 游 水 资源 的 开发 利 
用 和 浪费 ,来 水 逐渐 减少 , 台 特 玛 湖 湖面 水 位 不 断 
下 降 ,并 裸露 出 大 片 干 酒 湖 盆 区 ,湖面 面积 逐渐 变 
小 直至 干 润 ""。2000 年 以 来 , 随 着 塔里木 河流 域 综 
合 治 理工 程 建设 ,塔里木 河流 域 管理 局 先后 组 织 实 
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图 1 研究 区 地 理 位 置 及 区 域内 发 育 的 两 种 典型 风蚀 坑 
Fig. 1 Location of the study area and two typical blowouts 


developed in the region 


施 了 多 次 下 游 生 态 输 水 ,水 头 可 到 达 台 特 玛 湖 ,并 
可 形成 数 百 平方 千 米 的 湖面 ,此 时 的 湖面 基本 位 于 
218 国道 以 西 地 区 ,而 218 国道 以 东 的 原 干 润 湖 倪 因 
风沙 堆积 地 形 抬 高 ,生态 输 水 难以 到 达 ,依然 是 一 
片 风蚀 沙漠 化 快速 发 展 的 干旱 荒凉 之 地 。 
台 特 玛 湖 干 油 湖 盆 区 气候 极端 干旱 ,多 年 平均 
降水 量 仅 23.33 mm, 而 年 蒸发 量 达 2673.2 mm ,干燥 
引 数 为 63.0 ,平均 相对 湿度 为 41.1 匈 ;地 处 塔里木 盆 
地 东 灌 风 的 风口 地 带 , 感 行 NE ENE 风 , 风 力 强劲 ， 
多 年 平均 风速 为 2.16 ms', 多 年 平均 大 风 ( 风速 > 
17.2 m-s”) A ŽGK 28.2 d, 最 大 风速 17~25 m- s”, Wi 
时 最 大 风速 40 ms 7 ;而 河 湖 相 沉积 平原 的 沉积 
物 固 结 程度 较 差 ,结构 松散 而 干燥 ,地 表 植 被 稀少 ， 
JL BRERA ,部 分 地 段 发 育 了 厚度 不 同 的 盐 
壳 。 在 长 期 的 强风 沙 作 用 下 , 湖 盆 地 表 产 生 风 蚀 破 
口 ,下 伏 松 散 沙 层 出 露 ,水 平和 垂直 风蚀 速度 加 快 ， 
形成 大 小 不 一 的 风蚀 坑 以 及 多 种 形态 的 沙丘 。 该 
区 风蚀 坑 与 草原 地 区 风蚀 坑 不 同 呈 呈 ,处 于 平 沙 地 
中 , 坑 后 无 积 沙 区 ,是 由 于 风力 强劲 ,无 植被 阻挡 ， 
下 垫 面 较 草 原 地 区 粗糙 度 小 ,坑内 沙 粒 被 大 风 吹 
蚀 , 庄 挟 到 离 坑 较 远 地 区 。 经 卫星 遥感 影像 测量 发 
现 , 风 蚀 坑 类 型 主要 为 槽 形 和 碟 形 , 较 大 风蚀 坑 长 
轴 为 168m , 短 轴 为 70m, 较 小 风蚀 坑 长 轴 为 8m, 短 
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轴 为 Sm。 
2 数据 与 方法 


2.1 风 洞 试验 

实验 在 中 国 科学 院 新 疆 生 态 与 地 理 人 研究 所 黄 
索 湾 沙漠 研究 站 的 风 洞 实验 室 中 进行 。 该 风 洞 由 
中 国航 空 工业 空气 动力 研究 院 设计 和 建造 (图 2a)， 
属 直 流 吹 气 式 低速 风 洞 ,由 动力 段 、 整 流 段 .实验 段 
和 扩散 段 4 部 分 组 成 ,全 长 16.2 m, 其 中 实验 段 长 8 m 
(可 放置 模型 段 长 4m) ,截面 宽 1.3 m\、 高 1.0 m; 采 用 
侧 壁 扩 散 式 结构 , 侧 壁 扩散 角 0.2° ;实验 风速 0~20 
m's” 可 调 。 

实验 模型 用 石膏 制作 (图 2d~e) ,大 小 按照 实际 
比例 缩小 300 倍 ,具体 参数 如 图 24 中 所 示 。 实 验 时 
模型 表面 的 气流 速度 由 毕 托管 (高 度 分 别 为 1 cm 2 
em,3 em,7 cm、10 cm、15 cem、30 cem、50 cm) jill) gt JE 
力 计 算 而 得 (图 2c) ,气流 流向 由 粒子 图 像 测 速 仪 
(Particle image velocimetry , PIV, 型 号 Dynamic Stu- 
dio) 测 量 ( 图 2b)。PIV 可 以 跟踪 流 场 中 示 踪 粒子 
(直径 小 于 10 km, 甚至 am 级 ) 的 运动 ,并 用 数字 相 
机 同步 拍摄 流 场 照 片 ,系统 基于 前 后 两 帧 粒子 图 像 
来 计算 粒子 位 移 ,对 其 进行 互相 关 , 计 算得 到 流 场 
内 一 个 切面 的 速度 场 ,并 可 进一步 计算 得 到 切 应 变 
率 、 涡 度 等 各 种 运动 参数 。 


一 可 K = > | 
TEIN | 


将 风蚀 坑 模 型 置 于 风 洞 的 木质 底板 中 央 ,使 风 
蚀 坑 上 口水 平面 与 底板 齐 平 (图 2c)。 在 风蚀 坑 模 
型 中 轴线 上 布设 11 个 测 点 ,分别 在 6m.s ms '、 
12 mes "来 流风 速 下 测量 风蚀 坑 的 表面 风 场 。 基 于 
风速 数据 ,用 Surfer 10.0 绘 制 模型 表面 风 场 图 。 

用 玻璃 畦 将 激光 发 射 器 蛙 住 , 置 于 模型 的 正 前 
方 。 双 通道 脉冲 激光 采用 5 级 强度 ,间隔 160 mso 
高 速 粒子 摄像 机 置 于 模型 正 上 方 ,连续 拍摄 10 张 图 
片 , 以 记录 粒子 运动 方向 。 基 于 观测 数据 ,用 粒子 
流向 处 理 软件 (Dynamic Studio ) 绘 制 粒 子 流向 图 。 

在 33 cmx33 cm 沙盘 上 ,用 印 模 手 工 制作 大 小 
不 同 两 种 槽 状 风蚀 坑 沙 盘 模 型 ,大 模型 缩 尺 比 约 为 
1:320 ,小 模型 约 为 1:330 ,同时 设置 无 风蚀 坑 沙 盘 。 
分 别 在 7.5 ms-" 和 8.5 ms 风速 下 进行 净 风 吹 刨 实 
验 , 用 1/1000 kg 电子 秤 称 量 吹 蚀 前 后 沙盘 重量 , 计 
算出 风蚀 量 。 

2.2 数值 模拟 

使 用 模拟 软件 ANSYS 进行 流 场 模 拟 , 其 中 计算 
模型 采用 k-w 模 型 中 的 SST 模 型 (SST 模 型 集成 了 标 
准 k-s 模 型 和 k-o 模 型 的 优点 ,特别 适用 于 逆 压 梯度 
流动 和 分 离 流动 )。 本 文 模拟 的 风沙 流 马赫 数 均 小 
于 0.3 ,为 不 可 压缩 流动 ,计算 域 模型 出 口 边界 条 件 
为 压力 出 口 ,压力 差 为 0, 计 算 域 模型 人 口 边界 条 件 
设置 为 速度 入 口 , 和 人 口 为 虚 型 风速 廓 线 流 汪 ,大 小 


图 2 实验 设备 和 模型 


Fig. 2 Experimental equipment and model 
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为 33 mx33 mx11 m, PRAT A 22 m, 短 轴 
16 m, PRE 3.5 m; 碟 状 风 刨 坑 长 轴 短 轴 都 为 16 m, 深 
度 3.5m。 风 沙 流 中 沙 粒 粒 径 为 0.075~0.250 mm, FA] 
此 数值 模拟 中 沙 粒 粒 径 设 定 为 &=0.15 mm ,形状 为 
颗粒 , 沙 粒 密度 p.=2650 kgm”, if u.=0.0047 Pars, 
初始 沙 粒 体 积分 数 为 1% ,类 型 为 FLUID。 空 气 密度 
p=1.225 kgm”, ZAD J Bh FE = 1.78910 Pass, 
模拟 风速 选取 6 ms'。 计 算 网 格 采 用 三 角形 划分 ， 
网 格 数 量 分 别 为 365490、262050( 图 3), 质 量 良好 。 


(a) 槽 状 风蚀 坑 网 格 
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3.1 风 洞 实验 数据 分 析 

3.1.1 风蚀 坑 测 点 风速 廊 线 NERS RAI 
点 的 风速 廓 线 可 知 ,多 数 测 点 的 风速 廊 线 整体 呈 两 
段 式 ,在 3~7 cm 高 度 范围 内 其 至 出 现 风速 随 高 度 增 
加 而 降低 的 现象 (图 4); 而 风蚀 坑 外 的 测 点 15 风速 
RAG oT BO ,符合 正常 情况 ,说 明 测 点 15 气流 受 


(b) 碟 状 风蚀 坑 网 格 


图 3 网 格 划分 
Fig. 3 Mesh division 


(a) 各 测 点 风速 廓 线 
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图 4 风速 廓 线 图 
Fig. 4 Wind speed profiles 
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模型 影响 很 小 ,而 其 他 点 受 模型 影响 很 大 。 
3.1.2 风蚀 坑 的 纵 训 面 风 场 ”在 不 同 指示 风速 下 ， 
槽 状 坛 和 碟 状 坑 纵 剖面 风 场 图 可 以 看 出 ,模型 坑口 
以 上 主要 影响 高 度 在 20 cm 以 下 。 

在 6m:.s "指示 风速 下 (图 Sa) , 随 着 气流 靠近 风 
蚀 坑 内 背风 坡 处 ,两 种 模型 即 开始 呈现 不 同 的 风速 
变化 。 在 醒 状 模型 背风 坡 处 (X=11.4~15.9 cm), 出 
现 气流 加 速 现象 ; 而 在 碟 状 模型 背风 坡 (z15.3~ 
19.8 cm) 则 相反 ,出 现 气 流 减 速 现象 。 气 流 到 达 模 
型 中 部 时 ( 槽 状 了 15.9~22.2 cm HEAR X=19.8~24.8 
cm) , 槽 状 坑 等 值 线 上 凸 ,形成 减速 区 ,而 碟 状 坑 等 
值 线 下 四 ,形成 加 速 区 。 气 流 流 经 槽 状 风蚀 坑 尾 部 
迎风 坡 时 (X=22.2~33.2 em) ,等 值 线 下 四 ,形成 加 速 
区 ,而 流 经 尾部 迎风 坡 时 (X=24.8~29.3 em) ,7 cm 以 
上 等 值 线 上 凸 ,形成 减速 区 ,7 cm 以 下 等 值 线 下 四 ， 
仍 为 加 速 区 。 

在 9m:.s 指 示 风 速 下 (图 $b) , 槽 状 风蚀 坑 前 中 
部 (X=11.4~29.2 cm) 为 减速 区 ,尾部 迎风 坡 (X= 
29.2~33.2 cm) 为 加 速 区 , 较 6 ms' 指 示 风 速 相 比 ,加 
速 区 后 移 ; 碟 状 坑 风 速 呈 现 先 减 速 后 加 速 趋势 , 较 6 
ms"' 指 示 风 速 相 比 ,中 后 部 7 cm 以 上 变 为 完全 加 速 。 

在 12 ms! 指 示 风 速 下 (图 5c), 槽 状 坑 速 度 变 
化 趋势 与 9 ms 类 似 , 为 先 减 速 后 加 速 ,但 减速 区 
等 值 线 更 平滑 ,速度 递减 率 更 小 ; 碟 状 坑 速度 变化 
趋势 也 与 9 ms"' 类 似 ,但 等 值 线 最 凸 处 较 9 mesa 
移 1.5 cm, 且 减速 区 等 值 线 起 伏 大 ,速度 递减 更 快 。 
3.1.3 风蚀 坑内 的 气流 方向 图 6 为 图 像 粒 子 测速 
仪 获得 的 两 种 模型 表面 的 风沙 流 流向 图 。 由 图 可 


(a) 6 ms RUA 


(b) 9 ms"! 风 场 图 
25 


知 ,在 风蚀 坑 上 口水 平面 ,入 风 侧 处 气流 出 现 较 大 
变化 ,而 在 风蚀 坑 的 中 部 与 后 部 ,风沙 流 流向 与 初 

台 流向 一 致 。 在 入 风 侧 处 风沙 流 的 流向 受 模型 影 
响 而 发 生 偏转 (图 6 的 AC、E、G 区 ), 并 在 入 风 侧 处 
形成 风沙 流向 坑内 吸附 现象 ,也 有 风沙 流 从 坑 外 侧 
向 坑内 侧 聚 集 时 发 生 的 风沙 流 首 转 现象 (图 6 的 A、 
BCE、G 区 ), 形 成 了 以 模型 边缘 为 分 界线 的 对 向 风 
沙 流 (图 6 的 B 区 ), 而 碟 状 风蚀 坑 在 入 风 侧 (图 6 的 F 
区 ) 处 则 没有 风沙 流 逆 转 现象 。 
3.2 CFD 数值 模拟 数据 分 析 
3.2.1 模拟 真实 性 验证 在 以 上 风 洞 实验 与 数值 模 
拟 中 ,风蚀 坑 模 型 上 口水 平面 风速 具有 较 好 的 吻合 
性 (图 7) ,表明 k-o 模 型 中 的 SST 模 型 能 较 好 的 模拟 
风蚀 坑 的 动力 学 特征 ,这 解决 了 Smyth 等 ”所 使 用 K- 
8 模型 模拟 中 出 现 风 蚀 坑 背风 侧 不 准确 的 问题 。 
3.2.2 速度 通过 CFD 速度 模拟 发 现 , 在 两 种 风蚀 
坑 模 型 上 口 入 风 侧 背 风 坡 处 气流 首先 加 速 ,随后 减 
速 ,到 达 模 型 中 部 时 再 次 加 速 ,在 出 风 侧 处 形成 尾 
流 加 速 区 (图 8a、c) ,速度 变化 的 位 置 与 风 洞 实验 结 
果 一 致 。 两 种 风蚀 坑 模 型 差别 在 于 覃 状 模型 人 风 
侧 加 速 区 范围 大 于 碟 状 模型 。 从 模型 底面 可 以 看 
出 坑内 表面 流速 变化 情况 。 在 两 种 风蚀 坑 模 型 底 
部 ,气流 进入 风蚀 坑内 部 , 随 着 空间 增 大 气流 流 线 
发 生 分 离 ,在 入口 处 背风 坡 风 速 降 低 , 模 型 中 部 又 
开始 加 速 ,在 尾部 迎风 坡 处 大 幅度 减速 ,到达 模型 
外 形成 尾 流 加 速 区 (图 8b .d)。 

风蚀 坑 上 口水 平面 气流 与 底部 气流 大 小 不 同 ， 


(c) 12 ms-! 风 场 图 
25 


8 a 速度 /ms £ a 速度 /ms-! E E ems 
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图 5 模型 表面 风 场 剖面 图 


Fig.3 Profiles of wind field on the model surface 
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图 6 模型 表面 风沙 流 流 向 图 


Fig.6 Diagram of flow direction of wind and sand flow on the model surface 


6.5 F(a) 风 洞 实测 槽 状 坑 上 口水 平面 


图 7 实测 和 模拟 风蚀 坑 中 轴线 纵向 风速 变化 对 比 


Fig.7 Comparison between measured and simulated longitudinal wind speed changes in the central axis of blowout 


(a) 槽 状 坑 上 口水 平面 (b) 槽 状 坑 底 面 


(c) 碟 状 坑 上 口水 平面 (d) 碟 状 坑 底 面 


图 8 CFD 模 拟 速度 分 布 
Fig. 8 CFD simulated speed distribution 


上 口水 平面 气流 主要 受 空气 阻力 影响 ,而 底部 气流 
受 空气 阻力 和 地 面 摩擦 力 影响 ,因此 底面 流速 整体 
小 于 上 口 。 虽 然 两 种 模型 内 部 都 有 加 速 区 ,但 速度 
仍 小 于 初始 速度 (6 ms"'), 因 此 ,在 风蚀 坑 的 负 地 形 
作用 下 , 近 地 表 风 能 总 体 呈 被 削弱 状态 。 
气流 在 进入 风蚀 坑 后 ,发生 附 面 层 分 离 ,在 坑 
底 扩散 减速 ,并 在 模 状 坑 前 部 出 现 训 流 , 与 王 帅 等 ” 


和 孙策 等 ”野外 对 槽 状 坑 观测 结果 一 致 ,但 本 研究 
还 发 现 气 流 在 风蚀 坑 上 口 的 入 风 侧 处 气流 为 加 速 
状态 ,这 与 气流 集聚 辐 合 进 入 风蚀 坑 有 关 。 风 人 蚀 坑 
内 的 气流 变化 模式 不 仅 与 风蚀 坑 形 状 2 Ti] IN BE 
来 流速 度 有 关 ,还 受 风 的 入射 方向 影响 ”。 两 种 风 
蚀 坑 内 气流 变化 的 差异 主要 源 于 风蚀 坑 的 空间 太 
度 的 不 同 。 与 碟 状 坑 相 比 , 槽 状 坑 长 轴 更 长 ,气流 


202302.00029v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


1790 tid Ww 45% 


发 育 空间 相对 充足 ,速度 及 压力 变化 明显 。 

3.2.3 压力 ”风蚀 坑 模 型 表面 压力 分 布 能 反映 流速 
变化 和 气流 运动 方向 。 两 种 模型 上 口水 平面 (图 
9a.c) 都 在 入 风 侧 背风 坡 处 形成 低压 区 , 随 风沙 流 继 
续 深 入 ,压力 逐渐 增 大 ,在 模型 迎风 坡 处 达到 最 大 ， 
而 在 出 风 侧 处 因 风 沙 流 加 速 溢出 又 形成 低压 区 , 整 
体 中 部 压强 高 ,而 两 侧 压 强 低 ; 槽 状 模型 人 风 侧 低 
压 区 范围 大 于 碟 状 模型 ,而 碟 状 模型 后 部 低压 区 范 
围 较 大 (图 9c)。 模 型 底面 (图 9b、d) 压 力 变 化 趋势 
与 模型 上 口 不 同 ,压强 呈 前 低 后 高 分 布 模式 ,上 且 矶 
状 模型 底部 入 风 侧 低压 区 范围 更 小 ,而 后 部 没有 低 
压 区 。 在 底部 强烈 高 压 区 风沙 流 侵蚀 作用 下 , 坑 后 
缘 向 下 风向 扩展 守 ,这 与 风 洞 实验 观测 到 的 结果 一 


(a) 槽 状 坑 上 口水 平面 (b) HERS LR 


致 (图 10)。 此 外 图 10 中 发 现 有 微小 的 溯源 侵蚀 
现象 。 

3.2.4 风沙 流 流向 查 明 风蚀 坑内 风沙 流 流 向 有 助 
于 更 好 地 理解 风蚀 坑 的 形态 动力 学 过 程 。 图 11 为 
CFD 模拟 的 两 种 风蚀 坑内 气流 的 流向 ,可见 风 蚀 坑 
对 近 地 表 气流 产生 了 很 大 扰动 , 流 线 偏转 ,流向 趋 
于 复杂 化 , 坑 外 的 气流 受 局 地 气压 影响 ,被 吸附 到 
坑内 ,并 沿 侧 壁 行进 (图 1la.ce)。 在 风蚀 坑内 部 , 气 
流 运 动 更 加 复杂 ,在 入 风 侧 背风 坡 多 有 涡流 发 育 
(图 11b) ,而 在 模型 底部 中 心 线 处 ,气流 直线 前 进 ， 
随 着 向 两 侧 延展 ,气流 治 着 底面 坑 壁 呈 螺 旋 曲 线 前 
进 (图 11a) ,这 与 Hesp 研究 的 海岸 带 槽 状 风蚀 坑 气 
流 模 式 一 致 ” ,但 又 略 有 不 同 :气流 螺旋 前 进 现象 


(c) 碟 状 坑 上 口水 平面 


(d) 碟 状 坑 底 面 


图 9 _ CFD 模拟 压力 分 布 


Fig.9 CFD simulated pressure distribution 


(a) 吹 蚀 0 min 
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图 10 槽 形 模型 发 育 过 程 


Fig. 10 Slotted model development process 


(a) 槽 状 坑 上 口水 平面 (b) HEARSE 


x 
: 
1 


(c) 碟 状 坑 上 口水 平面 


图 11 CFD 模 拟 风沙 流 流向 


Fig. 11 CFD simulation of wind and sand flow direction 
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模型 人 风 侧 的 涡流 规模 尺度 小 于 档 状 模型 (图 
11e) , 矶 状 模型 侧 壁 风沙 流 沿 坑 壁 前 进 ,没有 发 育 涡 
流 ( 图 11d)。 风 蚀 坑 内 发 育 的 涡流 一 定 程度 可 以 解 
释 PIV 测 量 中 出 现 的 对 向 流 现象 。 


4 讨论 


4.1 风沙 流 吸 附 机 制 

以 上 模拟 发 现 , 近 地 表 和 气流 对 两 种 风蚀 坑 地 形 
表现 出 极其 复杂 的 响应 ,尤其 发 现在 权 状 坑 和 碟 状 
坑 的 坑口 入 口 处 都 有 加 速 现象 。 这 与 周围 气流 向 
负 地 形 辐 合 汇 集 有 关 , 而 气流 汇聚 则 与 气流 自身 重 
力 和 坑口 局 部 的 低压 区 有 关 。 气 流 流 问 模拟 也 发 
现 风 刨 坑 具 有 吸附 效应 , 流 经 坑 外 一 定 范围 的 气流 
被 吸附 到 坑内 (图 11、 图 12) ,坑口 及 其 上 方 一 定 高 
度 范 围 的 气流 都 向 坑内 汇聚 (图 12)。 测 量 发 现 , 在 
6 m-s KUET ,两 种 风蚀 坑 两 侧 吸 附 范围 大 约 为 短 
半 轴 的 0.25 倍 ,垂直 吸附 范围 大 约 为 深度 的 1.56 
倍 , 且 多 发 生 于 风蚀 坑 的 前 中 部 。 由 伯 努 利 效 应 


(a) 槽 状 坑 气 流 吸附 
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(不 可 压缩 流体 流速 越 快 ,流体 产生 的 压力 就 越 小 ) 
可 知 ,坑口 处 流速 增 大 导致 压力 变 小 ,因而 使 周围 
气流 在 压力 差 作用 下 被 吸 和 人 坑内。 根据 流体 力学 
的 能 量 守 恒定 律 , 由 于 风蚀 坑 的 吸附 机 制 ,两 侧 和 
上 方 的 气流 携带 能 量 被 聚集 到 坑内 ,使 坑内 气流 能 
量 增加 ,同时 也 将 风沙 流 带 入 坑内 ,而 风蚀 坑 出 风 
侧 处 气流 能 量 较 集 中 释放 (图 8) ,加 剧 了 风蚀 坑 出 
风 侧 的 侵蚀 ,使 风蚀 坑 逐 渐 发 育 。 孙 画 等 汪 在 槽 状 
坑 气 流 观测 中 也 观测 到 气流 被 吸入 坑内 的 现象 ,但 
没有 量化 被 吸入 气流 的 范围 ,也 没有 从 能 量 汇聚 角 
度 分 析 对 风蚀 坑 发 育 的 意义 。 本 研究 较 全 面 获取 
了 风 刨 坑 的 流 场 信息 ,证实 了 涡流 的 存在 ,并 阐明 
了 风蚀 坑 的 风沙 流 吸 附 机 制 , 有 助 于 深入 理解 风蚀 
坑 的 风沙 动力 过 程 。 

风蚀 坑 的 吸附 机 制 汇聚 能 量 加 剧 对 地 表 的 侵 
蚀 。 对 不 同 沙 盘 采 用 不 同 风速 吹 蚀 12 min 后 发 现 
( 表 1) ,在 大 风 条 件 (8.5 mes") FD) RE Ry 
盘 和 大 尺度 风蚀 坑 沙 盘 的 风蚀 量 百分比 分 别 比 平 
坦 沙 地 沙盘 的 大 3.14% .6.96% ;在 小 风 条 件 下 DIR 
度 风 蚀 坑 沙盘 和 大 尺度 风蚀 坑 沙 盘 风蚀 量 百分比 


(b) 碟 状 坑 气流 吸附 


-52.2 


12 风蚀 坑 流 场 显示 的 气流 吸附 效应 
Fig. 12 Adsorption effect of airflow in blowout flow field 


表 1 风蚀 坑 对 地 表 蚀 积 的 影响 


Tab.1 Influence of blowout on surface erosion and deposition 


影响 因素 7.5 mes -风速 8.5 m- s” KUE 
平坦 沙 地 沙盘 小 尺度 风蚀 坑 沙盘 大 尺度 风蚀 坑 沙盘 ”平坦 沙 地 沙盘 小 尺度 风蚀 坑 沙 盘 大 尺度 风蚀 坑 沙盘 
质量 /kg 吹 蚀 0 min 6.526 6.526 6.526 6.526 6.526 6.526 
吹 蚀 12 min 6.337 6.321 6.309 5.506 5.474 5.435 
风蚀 量 /kg 0.189 0.205 0.217 1.020 1.052 1.091 


风蚀 量 百分比 /9 100.00 108.47 114.81 100.00 103.14 106.96 
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分 别 比 平 坦 沙 地 沙盘 的 大 8.47% 14.81%。 可 见 ,在 
相同 风速 下 有 风蚀 坑 沙 盘 的 风蚀 量 大 于 平坦 沙 地 
沙盘 , 且 风 蚀 坑 尺度 越 大 风蚀 量 越 大 ;相同 模型 条 
件 下 风速 越 大 ,风蚀 量 越 大 。 因 此 ,侵蚀 程度 与 风 
速成 正比 ,风蚀 坑 加 剧 了 地 表 侵 蚀 , 且 侵蚀 程度 与 
风蚀 坑 尺 度 有 关 , 风 蚀 坑 尺 度 越 大 侵蚀 越 强 。 
4.2 风蚀 坑 发 育 的 关键 机 制 

前 人 研究 都 已 证 实 风 蚀 坑 的 起 始 和 发 育 源 于 
外 部 因素 (风沙 流 呈 沙丘 后 ) 和 内 部 因素 ( 龟 裂 纹 ) 
的 耦合 "。 当 地 表 产 生 风 蚀 破 口 时 ,微小 的 负 地 形 
对 局 部 气流 运动 产生 了 阻 清 ,气流 经 历 扩散 减速 和 
溢出 加 速 过 程 , 破 口内 的 沉积 物 随 之 被 侵蚀 搬运 ， 
破 口 逐 渐 扩 大 和 加 深 ; 随 之 对 气流 阻 滞 作 用 也 加 
强 , 下 伏 沉 积 物 被 气流 侵蚀 吹 走 ,形成 表面 光滑 的 
流线型 外 形 ,形成 雏形 坑 。 不 同 风向 的 风沙 流 对 纵 
形 坑 可 从 不 同方 向 进行 侵蚀 和 形态 塑造 ,并 最 终 形 
成 与 区 域 风 况 相 适应 的 外 形 ; 委 形 坑 规模 逐渐 扩 
大 ,风蚀 坑内 气流 的 阶段 性 变化 更 加 强烈 , 潜 蚀 强 
烈 发 育 ,使 下 伏 更 多 沉积 层 出 露 ,形成 质地 软 便 相 
间 的 侧 壁 ,并 产生 强烈 的 差异 性 侧 刨 ””, 在 重力 作 
用 下 发 生 块 体 册 塌 , 侧 蚀 快 速 发 展 ,逐渐 形成 碟 状 
坑 和 槽 状 坑 ; 风 蚀 坑 逐步 扩大 ,相互 联通 ,形成 风蚀 
谷 和 风蚀 残 墩 部。 在 台 特 玛 湖 干 泗 湖 盆 区 , 覃 状 坑 
和 蝶 状 坑 形 态 差 异 与 成 行 风 方 向 、 发 育 时 间 、 下 热 
面 植被 状况 及 微 地 貌 起 伏 状况 等 有 关 。 单 一 风 况 
下 ,发 育 时 间 较 长 ,成 为 槽 状 坑 ; 相 对 复杂 风 况 , 且 
发 育 时 间 较 短 ,成 为 碟 状 坑 。 权 状 风蚀 坑 常 有 陡峭 
的 侧 壁 ,而 在 坑 底 沿 侧 壁 有 少量 崩塌 的 块 体 ™, 这 
就 是 气流 对 软 人 硬 相 间 侧 壁 沉积 物 侧 向 差异 性 掏 蚀 
的 表现 。 侧 壁 掏 蚀 的 气流 有 沿 纵 轴 进 入 坑内 的 侧 
向 绕 流 和 螺旋 涡流 ,也 有 坑 体 周边 吸附 进入 的 
涡流 。 


5 结论 


li 


Sa pst St ss WS Pe Ss X AI ES HH AUT AS, tH ee DAL 
foe ee AS BO AFR RWT 
TTAB 2 K EA BON i >, RA EY USI ASB 
别 于 草原 和 海岸 区 , 仅 有 坑 体 ,而 没有 坑 后 积 沙 区 ， 
形态 以 权 状 坑 和 碟 状 坑 为 主 。 风 洞 观测 和 CFD 模 
拟 显示 : 

(1) 两 种 风蚀 坑 对 近 地 表 气流 造成 强烈 扰动 ， 


权 状 风蚀 坑 和 碟 状 风蚀 坑内 气流 速度 和 压强 变化 
模式 基本 相似 ,坑口 气流 速度 变化 呈 加 速 一 减速 一 
加 速 , 压 强 旦 低压 一 高 压 一 低压 ,在 出 风 侧 处 形成 
加 速 尾 流 恢复 区 ,而 底面 呈 减 速 一 加 速 一 减速 一 加 
速 , 压 强 呈 低压 一 高 压 , 在 出 风 侧 处 形成 高 压强 气 
流溢 出 区 ,两 种 模型 气流 差异 在 于 个 别 变 化 区 范围 
的 大 小 。 

(2) 研究 发 现 ,由 于 负 地 形 对 气流 的 扰动 ,在 入 
风 侧 处 发 育 涡流 ,在 覃 状 坑 侧 壁 有 螺旋 环流 ,内 部 
有 弧 形 前 进 流 。 

(3) 两 种 风蚀 坑 均 具有 一 定 强度 风沙 流 吸附 歼 
应 ,将 水 平和 上 层 气 流 的 能 量 和 风沙 流 聚 集 到 坑 
内 ,加剧 了 坑 体 的 风蚀 速度 。 风 蚀 坑 吸 附 效应 是 风 
蚀 坑 发 育 的 关键 作用 机 制 。 从 风 洞 实验 研究 也 发 
现 ,风蚀 坑 加 剧 了 地 表 侵 蚀 ,侵蚀 程度 与 风蚀 坑 尺 
FEAT ,风蚀 坑 尺度 越 大 对 地 表 侵 蚀 越 强 。 
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Three dimensional flow field characteristics of two typical blowout in the dry 
lake basin of Taitema Lake and their influence on wind erosion 
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Abstract: Blowouts in the dry lake basin of Taitema Lake, south Xinjiang, China are primarily slotted pits and 
saucer- shaped pits, and the surface wind erosion is not only influenced by the incoming wind velocity but also 
strongly influenced by the three-dimensional (3D) flow field in these blowouts. Based on wind tunnel experiments and 
computational fluid dynamics numerical simulations, the 3D flow field of blowouts is investigated. The following 
results are found: (1) From the air inlet to the air outlet, the wind speed at the bottom of the grooved and dished 
blowouts exhibits a change mode of deceleration-acceleration-deceleration- acceleration. The wind speed in the 
horizontal plane, where the upper mouth is located, exhibits a change mode of acceleration-deceleration- acceleration. 
The air pressure at the upper mouth also changes spatially accordingly. The overall change mode of the grooved and 
dished blowouts is low pressure-high pressure-low pressure. The change of trough pit is more obvious. (2) Both 
blowouts have a certain adsorption effect on wind and sand flow, causing wind and sand flow into the pits, 
accumulating energy, wind, and sand flow and intensifying wind erosion development. (3) Blowouts intensify surface 
erosion, and the degree of erosion is related to the scale of the blowout; the larger the scale of the blowout, the stronger 
the surface erosion. This study visualizes the 3D velocity, pressure, and wind-sand flow distribution patterns in trough- 
shaped and disc-shaped blowouts and proposes a wind-sand flow adsorption mechanism in these blowouts, the results 
of which can provide theoretical guidance for inland dry lake basin desertification control. 


Key words: blowout; wind tunnel experiments; CFD numerical simulation; wind and sand flow 


